NPS ARCHIVE 
1963 

HEERWAGEN, D. 



INTERSTELLAR MATERIAL IN SLOSULAR CLUSTERS 
DAVID DUPUY HEERWAGEN 




LIBRARY 

US. RXval postgraduate school 

MONTEREY, CALIFORNIA 



DUDLEY KNOX LIBRARY 

S^v 1 ?? A S JATE school 

MONTEREY CA 93943-5101 



INTERSTELLAR MATERIAL 



IN 

GLOBULAR CLUSTERS 



by 

David Dupuy Heerwagen 
LCDR,, U.S. Navy* 



Thesis submitted to the Faculty of the Graduate School of 
the University of Maryland in partial fulfillment 
of the requirements for the degree of 
Master of Science 
I960 



*The author attended the University of Maryland under the sponsorship of the 
Uo So Naval Postgraduate School,, The assertions contained herein are those 
of the author and are not to he construed as official or reflecting the views 
of the Navy Department or the naval service at large „ 



he bis 



'a/ , 0 ' 




library 

U.S. NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 



APPROVAL SHEET 



Title of Thesis: Interstellar Material in Globular Clusters 



Name of Candidate: David Dupuy Heerwagen 

Master of Science,, 1963 
/ / . 






ABSTRACT 



Title of Thesis: Interstellar Material in Globular Clusters 

David Dupuy Heerwagen* Master of Science* 1965 

Thesis directed by: Uco van Wijk 

Assistant Professor 

The chances of survival of interstellar material in globular clusters is 
investigated. The observational evaluations of Roberts (i960) and Idlis and 
Nickol'skii (1959) used. The effects of ultraviolet radiation due to 
hot* blue stars in the center of the cluster and the effects of the radiation 
field in general are discussed in Chapter III. The reaction due to collisions 
of* globular cluster clouds with material in the galactic disk is evaluated 

both from momentum and energy influx considerations* and the probability of 
ejection of clouds from, the cluster is estimated in Chapter IV, It is sug- 
gested that the most serious uncertainty is introduced by the incompleteness 
factor for hot 0 stars in globular clusters and further observations of these 



not stars are desired* 
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CHAPTER I 



INTRODUCTION 

Two recent investigations have stirred again the controversy concerning 
the possible presence of obscuring matter in globular clusters. First* Idlis 
and Nikol'skii (1959) reported sighting tiny globule-like objects in the central 
region of the globular cluster closest to the sun* Messier 4. Roberts (i960) 
investigated obscuring* cloud -like objects which were" discovered in 12 out of 
32 high altitude globular clusters. 

The arguments presented by these investigators that the obscuring material 
is an integral part of the globular clusters seem convincing. They both showed 
that the probability of such vacant or obscured regions appearing as a result 
of statistical fluctuations in stellar distribution of the cluster is extremely 
remote* if not impossible. Roberts also demonstrated that for such obscura- 
ion to be caused by intervening halo clouds would require a space density of 
uch clouds much higher than in the solar neighborhood - a result again extremely 
unlikely. It was also noted that clouds or obscurations do not appear against 
the background of the elliptical galaxies. 

In this investigation then^ we start with the assumption that these objects 
are indeed obscuring medium in the clusters* to be analyzed concerning both 
their reaction with the galactic disk and the possible environment within the 
cluster which affects their formation and behavior. From the investigation 
of Roberts* we sense a possibility that* since not all the clusters show obscur- 
ing matter* there is a selection effect operating to make conditions for forma- 
tion of such objects and their retention and persistence within the clusters 
more favorable in certain cases than in others. This approach then suggests the 
possibility that there may be small but significant differences between individual 
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clusters. However., the globular clusters seem to he a remarkably homogeneous 
group of objects with respect to certain properties* Von Hoerner’s (195?) study 
of globular clusters showed that their standard deviation in luminosity is less 
than 1 absolute magnitude., that their deviation of diameters is approximately 
2 6 jo of the mean diameter * and that their internal structures are very similar* 
This homogeneity suggests that we adopt a standard cluster model as the frame- 
work within which to make our order of magnitude calculations concerning the 
interstellar material* Because of the great amount of observation and study 
devoted to the cluster Messier 3^ we have generally adopted this object as our 
model* 

Of the many basic assumptions that must be made in an investigation of 
this sort., perhaps the most difficult is that concerning the physical state 
and concentration of matter in the galactic disk^ lying outside the discrete 
clouds and galactic spiral arms. Results from, the study of optical,, inter - 
stellar absorption lines,, emission nebulae., absorption and obscuration of 
starlight «, and from radio astronomy work in both the continuum and with the 
21 -cm line have given a fairly consistent picture of the galactic structure 
in the dense regions* However., a search of the literature reveals that con- 
fusion and controversy presently exist on the state and density of matter 
between the arms. 

In his study of absorption lines in the optical region,, Stromgren (1948) 
implies that the non-cloud region is permeated with a HI medium of approxi- 
mately 0.1 H atoms cm 0 Spitzer (19^-9) suggests at the first symposium on 
eosmical gas dynamics that the interstellar medium is in pressure equilibrium 
except for the relatively small volume of expanding H II regions near hot stars. 
His model would require the diffuse intercloud medium to have a kinetic tempera- 
ture of 10 j 000 °K<, signifying that it is primarily ionized or HII. A summary 
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of his model shows s 
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However, later Spitzer (195*0 states that the "atoms in this tenuous intercloud 

medium could be mostly neutral and still at a relatively high temperature", Dufay (1955) 

-2 -3 

says that regions of diffuse hydrogen of 10 cm concentration could be HI regions 
of temperature 1000 K or less., because of the long times (l(r years) required to 
reach equilibrium at this density,, A gas at this density cools at a very slow 
rate. However, on page 123 he remarks that much of the intereloud medium could 
be Hilo Finally,, the question of whether the medium is predominantly HI or HI I 
is discussed at length in the Third Symposium on Cosmical Gas Dynamics held in 
1958- Spitzer (1958) concludes that "in any case,, this intercloud medium is 
probably at a high temperature, so thau its pressure is equal to that in the 
censer, cooler HI regions. Whether this rarified gas between clouds,, and not 
near an early type star,, is HI or HII is an open question". 

In all probability., the intercloud hydrogen of the galactic disk is not uni- 
form in concentration, temperature, or degree of ionization. However, an upper 
limit to its concentration can be set reasonably as 0,1 cm or one order of magni- 
tude lower than the concentration of 1 H atom cm assigned to the spiral arms by 
radio 21-cm line and optical, absorption line studies. Concentrations as low as 
0,01 cm can also be considered reasonable. 

In this investigation, we are concerned with the gross dynamics of the encoun- 
ter process between a dense cloud or globule of the globular cluster and a typical 
cloud or the intercloud medium of the galactic disk. To simplify our calculations 
and to enable us to obtain both order of magnitude answers and limiting parameters, 
we have assumed the intercloud medium to be neutral. In fact, by using the 
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macroscopic approach of transfer of momentum and kinetic energy between large 
regions of gas and ignoring the individual particle or microscopic behavior, 
we reason that the question of a neutral or an ionized condition for the inter- 
cloud medium, interacting with the neutral clouds of the cluster, is not a 
sensitive input to the final estimates. Bates and Griff ing (1953) have calcu- 
lated collisional ionization cross sections for the neutral -neutral and neutral- 
ionized interaction of hydrogen. They find, that at the lower energies of approx- 
imately the range of our problem, the ionization cross section of the two reac- 
tions are approximately the same. 

Little of theory or experiment, applying to astrophysical problems of 
these- densities, has been done for the case of plasma interactions, Alfven 
(i960) has considered the case of neutral -ionized collision. He analyzes the 
case where a non-ionized gas collides with a plasma where the kinetic energy 
per particle of relative motion between the particles of the two media exceeds 
the ionization energy per particle of the neutral gas. He concludes that the 
neutral gas will be ionized and trapped in the magnetic field and the density 
of the plasma will increase. However, implicit in his treatment, and also in 
the two papers by Kahn (1957? 1958) concerning collision of two ionized gases, 
is the assumption that the two gases have comparable densities. We will find 
that the case which ultimately interests us in this investigation is the colli- 
sion of a dense, neutral gas of relatively small volume with a highly tenuous 
gas of much larger extent, 

Kahn finds when two rarified masses of non-ionized gas collide at high 
enough velocity that both gases will be ionized after a relatively shallow 
penetration, and that the counter streaming electrons will be stopped in a very 
short distance. This results in a stationary, oscillating electron gas through 
which protons are streaming in opposite directions. This so-called collective 
instability of the electrons and not the individual collisions between particles 
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is .responsible generally for stopping the counter -streaming gases * 

However, Kahn does assume equal density of particles in each stream,, For 
the case where a wide difference exists in the concentrations of the two counter - 
streams, a momentum exchange treatment is appropriate to determine the result- 
ing motion of the arrested, now unified conglomerate. 

So we will restrict our attention in this problem to the question of what 
happens to the cloud of a globular cluster with respect to the cluster 1 s 
standard of rest when it impacts with the medium of the galactic disk. We will 
neglect considerations of individual particle reaction. 
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CHAPTER II 



DISCUSSION OF OBSERVATIONAL DATA AND ASSUMPTIONS 

Idlis and Nikol'skii (1959) detected small obscurations in the central 
regions of M4. These globules appear twice as numerous photographed in the 
blue as they do in the red,, thus following a rough l/\ absorption law and indi- 
cating the presence of dust particles characteristic of the plane of the Galaxy. 

They estimated the following average properties for the globules: radius = 0.02 pc., 

optical thickness (x) = 3, mass of dust = 0.001 solar mass (<2>), total mass = 0.10,, 

5 -5 

concentration of hydrogen atoms = 10 cm , total number of clouds = 140. They 
report these globules as concentrated toward the center of the cluster with maxi- 
mum concentration at about 0.6 pc distance from the center. 

On the other hand, Roberts (i960) has reported the detection of larger., 
more massive objects in 12 out of 32 high latitude globular clusters examined. 

A typical cloud is that reported in M 13 as having a projected image of 1 by 3 pc 
and a total photographic absorption of 5 magnitudes. Roberts estimates the mass 
of dust in this particular cloud as 2.6 O, Assuming the estimated 1$ ratio of 
dust to gas in the Galaxy plane to hold for the globular clusters., he finds the 
total mass to be 260 O. Actually, for simplicity of calculations, this is as- 
sumed to be the weight of the atomic hydrogen atoms. This would amount to a 
hydrogen atom concentration of the order of lCr cm . He estimates that a typical 
cluster which shows these clouds would have about 1000 O tied up in such objects. 
Roberts does not discuss the distribution or location of his clouds, other than 
to show by photograph that they generally appear at the cluster edge where the 
apparent density or luminosity appear to drop sharply. However, from his discus- 
sion of the dimensions of a cloud shown in a photograph of M 13, one can deduce 
that the cloud lies a minimum distance of about 8 to 10 pc from the cluster center. 
The procedure used to determine the 5 magnitude absorption is not discussed by 
Roberts and is probably a very rough estimate. 
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It is clear that the foregoing two investigators are concerned each with 
an entirely different object,, differing by about two orders of magnitude in 
both density and radius. Accordingly, any discussion of the possible origins, 
lifetimes, or evolutions of interstellar material in globular clusters should 
concern itself and attempt to account for both classes of objects. Throughout 
this paper, those objects reported by Idlis and Nikol'skii will be called 
globules while those of Roberts will be termed clouds. 

In making our analysis, we will adopt a detailed model of a typical globu- 
lar cluster. The cluster. Messier 3, has been observed and studied in great- 
detail by Sandage (l95*b 195?) and others. Accordingly, we use this observa- 
tional data as a basis for cur study. Also, we adopt Oort and van Herk's (1959) 
dynamical model (d) of M3 since it corresponds so closely to its observed structure. 

Interstellar material could come from any or all of three sources: (l) resi- 

dual material left over after original condensation of the cluster, (2) material 
that has been ejected from evolving cluster stars, (3) material that has been 
picked up from the halo and galactic disk by the accretion of clouds within 
the cluster. As we see later, sources (l) and (3) contribute a negligible 
amount compared to (2). We adopt the results of Sandage for the evolution of 
M3 which show that approximately 100,000 O have been converted into interstellar 
material over about 10^ years lifetime. We also assume the material to have 
the cosmic proportion of elements, as adopted for interstellar material in the 
galactic disk. 

Orbits of Globular Clusters - Most of the rough orbit calculations of globular 

clusters are based on the radial velocity measurements made by Mayall (19^-6) or 

by Kinman (1959)- In this investigation we will use the more recent data furnished 

by Kinman. He finds that the clusters follow highly eccentric orbits, the average 

eccentricity being 0.8. It is interesting to note that a graph in a recent paper 

by Eggen, Lynden-Bell, and Sandage (1962) also shows an average eccentricity of 

about 0.8 for high velocity stars of the subdwarf variety, which from ultraviolet 
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excess measurements , appear to correspond to the main sequence of globular 
clusters. 

Matsunami and collaborators (1959) have computed galactocentric distances 
for 9k of the 118 known clusters , using Lohmann's heliocentric distances. We find 
the average galactocentric distance to be 8.5 kiloparsecs (kpc) for their 9k 
clusters. Since we know that for highly eccentric orbits, an object is most 
apt to be found near apogalacticon, we make a rough estimate that the average 
apogalacticon distance (Q) is 10 kpc. Using the formula Q = a(l+e), we find 
the semi -major axis to be about 5-5 kpc. Using the formula for the latus rectum 

p 

p = a(l-e ), we find p = 2 kpc. 

Having found the parameters for the average orbit, we use them to obtain 
the average velocity (relative to the galactic frame of referneee) of a cluster 
as it traverses the galactic disk. For this calculation we assume that the 
average orbital axis makes an angle of V? degrees with the galactic disk, and 

we use the polar equation for an ellipse, r = - (where ^ = anomaly), 

solving for r at the two points of intersection. Upon finding r = k„5 kpc and 
1.3 kpc, we select r = 3 kpc as a typical intersection point on the galactic 
disk. 

We use the Schmidt model galaxy (1956) for the distribution of mass in the 
Galaxy, We find the following values from his graph for potential (in km /sec ) 

p p 

in the (w, z) plane; at Q = 10 kpc, = 30,000 km / sec ; at r = 3 kpc, 

2 2 

Sy = 80,000 km / sec . Adopting the value Vq = 60 km/sec at 10 kpc and using 
the formula = ^V r - we find V = 330 km/sec. 

We have made the above approximations for a "mean” orbit for a globular 
cluster, in order to arrive at some typical value for the velocity of a cluster 
(relative to the galactic center) as it traverses the disk. In reality the 
orbits are not closed as we will see in Chapter V, However, for the high e orbits 
with which we are concerned, the foregoing simplifications will give an adequate 
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estimate. We intend to use this typical value in our analysis of the probable 
reaction between clouds of the cluster and the interstellar medium of the disk 
during transit. We must recognize that any such analysis is by necessity, 
order of magnitude accuracy. We must realize that for each individual transit, 
since the disk clouds rotate with the circular velocity of the disk, the rela- 
tive velocity for collision can vary from about zero to approximately 600 km/sec. 

For purposes of our analysis,, we will select 300 km/sec as the mean relative 

8 

velocity. We use 2 x 10 years as the period for the mep,n orbit. 

Stability of Clouds and Globules - We have adopted Roberts's cloud as being a 
spherical object of 260 solar masses and 2 pc diameter. It is interesting to 
consider the gravitational stability of such an object. For this purpose we 
use Jean's gravitational instability criterion (1928), and find a rough limit 
that a self gravitating objecu must satisfy in order to be gravitationally 
stable and not break up under the result of its own thermal motions. In apply- 
ing any such criterion, we assume the cloud to be isolated in space. If there 
were a surrounding medium an a finite pressure, then the minimum temperature 
for stability would be larger. The most convenient form of this criterion, 
which can also be derived from the virial theorem, is: N>Ncr = 2.5x10^. T^/m^ 

This signifies that a body of mass M (in solar masses) at temperature T must 

-3 

have a density of more than cm ^ of hydrogen atoms to be stable. For the 
2o0 O cloud, we find the density of hydrogen atoms to be about 2700 cm . For 
tnis cloud to be gravitationally stable, its temperature would have to be below 
^0° to 50 °K. Such a low temperature is certainly feasible in an HI region 
( neutral hydrogen) particularly when appreciable obscuring dust is present, 
and there is a strong possibility that any such object is stable. On the other 
hand, the 0.1 O globules fall short of the critical density. However, Idlis 
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has shown that these objects should have a lifetime of approximately 10^ years 
when considering their thermal diffusion^ so the globules are probably long- 
lived objects in their present state* 
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CHAPTER III 



ENVIRONMENT WITHIN THE GLOBULAR CLUSTER 

If, in fact, we are convinced that the clouds and globules do exist in 
certain globular clusters, then the next matter of interest is an analysis of 
the environment within these clusters which encourages the formation and 
existence of such objects,, Of principle interest is an estimate of the amount 

o 

radiation from the stars beyond the Lyman limit (912 A), which would ionize the 
hydrogen. Again, we use the observational data available for M3. 

Star counts by Sandage for M3 ( 195^0 have shown that the group of stars brighter 
than absolute visual magnitude +3.5 ( 1.25 ©) is slightly concentrated toward the 
center of the cluster. It is in this group of more massive stars above the main 
sequence that we expect to find those actively evolving toward the white dwarf 
stage by ejecting material. We might, therefore, expect to find this matter ini- 
tially with the same distribution as these bright stars. Generally, this would 
mean a spherical distribution, somewhat concentrated toward the center. More 
specifically, we can give Roberts's estimate that a maximum of about 1000 solar- 
masses of interstellar material is available between passages through the galactic 
disk. He makes this estimate for a recent epoch by first finding the number of 
cluster stars (AN) lying above the main sequence for M3, using the luminosity 
function of Sandage (1957)“ He also computes the time (At) necessary for the star 
to evolve through this region of the color -magnitude diagram. Then after finding 
4*0, the average mass loss suffered by a star during evolution, (using 1.2 O for 
initial mass and 0.6 O for final white dwarf mass) he estimates the mass loss of 
the stars of the cluster for a half period of its galactic orbit (-|t) by the 
product; AN . ( -'■'^7 ) • " 3 ^. Since more than half of these bright stars lie 
within 8 pc of the center of our M3 model, we assume as a first approximation that 
we have 700 O of diffuse material spread evenly within this radius. 

This approximation allows us to make use of the concept of "Stromgren Spheres" 
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from Stromgren (1939) to show the extent of ionized hydrogen around stars of types 

05 to AO, assuming that the surrounding material is distributed uniformly and of 

-3 

density 1 H atom cm . These tabulated radii may then be modified for regions 

-2/3 

of different concentration by multiplying by N ' , 

To make this estimate for our cluster models we use the computations of 
Minkowski and Osterbrock (1939) for M3» They found that about 15 hot blue 
stars, of effective spectral type 05 and of considerably fainter absolute magni- 
tude (M = +2.3) than ordinary Population I stars of the same effective temper- 

* 

attire , would be responsible for most of the ionization. The uncertainty of the 
data is illustrated by the fact that only one of these 05 stars has been sighted 
and the rest are inferred from an incompleteness factor computed by Sandage (1954). 
These 15 stars are distributed within 28 parsecs of the center of the cluster. 

Since we are concerned in our calculation with the volume within 8 pc of the 
center, we correct for the lower number of stars by noting that for the Oort- 
van Herk model (d), the volume within 8 parsecs contains about 0.6 of the mass 
of the stars brighter than M = +3° 5 contained within 28 pc. So we correct the 



Minkowski -Osterbrock equation for the final total ionized volume, V = J-°9 X !0 (pc^) 

6 N 2 

where N number density, cm , of H atoms, and use V = — — ( p ) . Stromgren 

IT 

(1948) has given a correction to his tabulated radii of ionized spheres for 
the dilution effect for stars that lie outside the interstellar medium. Using 
this correction, we find we can neglect the contribution of those blue stars 
outside of 8 pc for our order of magnitude estimate. 

Now, using a total average mass of hydrogen atoms in diffuse matter within 
8 pc as 700 O, we find the average concentration of H atoms is N = 15 cm“^. At 
this density, we find that the whole volume within 8 pc would be ionized. This 
result is for a uniformly distributed medium, and does not take into account that 
any material would be concentrated toward the center as indicated by the potential 
field in Figure 1. Also there could and probably would be turbulence and density 
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fluctuations. Nonetheless , this particular environment does not appear conducive 
to the formation of HI regions., particularly of the size of 260 O. Ionized 
hydrogen has a kinetic temperature of approximately 10,000 °k, which gives a 
random thermal velocity of about 1 6 km/sec that would cause it to escape from 
the cluster quickly. 

In order to form a general impression of the effect of more interstellar 
material within a globular cluster, we will hypothesize that 3,000 O of mater- 
ial is available., which would be the case if the material remains ■ undisturbed in 
the cluster foi* & time interval • equal- to t'iiree mean orbital periods'. 

About 2,000 O of this would be within 8 pc. We will again assume it to be 

spread uniformly. Using the relation V = L^i^L (pc^) with N now equal 42 cm"^ 

JT , 

we find that less than one-third of the volume within 8 pc., or 630 pc" , is now 
ionized. We could visualize,, in this particular case., regions of ionized mater- 
ial or emission nebulae surrounding each hot,, blue star,, separated by larger 
volumes of HI gas. For the case between the two extremes just discussed,, we 
might picture instead a region of tenuous, ionized gas escaping around floating, 
smaller volumes of HI. 

The 2., 000 O of interstellar hydrogen would be just ionized, if distri- 
buted uniformly within a volume of 11 pc radius. It is interesting for this 
case to estimate the time interval for ionizing the entire 2,000 O, assuming 
that all quanta beyond the Lyman limit (912 A) are fully absorbed by the hydrogen 
atoms. To make this approximation, we must first calculate the radius of the 
sub-luminous 0-5 type star, which is the principal ionizing agent. We use a 
formula from Russell, Dugan, and Stewart (1938) My = ^2^22 _ 5 (qq R-0.08. 

From Strcmgren, we adopt an effective temperature for the 05 star of J^pOOO °K. 
Solving for R, we get 0.4 solar radii. 
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The energy density at the surface of a star of effective temperature, T, 



is given by the Planck law asj 

BIT h V 3 
Li V 



cP (e^T-i) 

(V ) 

To find the energy density in the frequency range beyond the Lyman limit' 
we integrate the above expressions 

u 



= (e h ^ r - lj'V 3 dV 
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At these high frequencies we make use of the Wien approximations f — Q 
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Integrating, we get; 
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Now, to find the intensity at the star's surface, we use the relation: 
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This is the luminosity beyond the Lyman limit for one 05 star of M y = +2.3» 

We can check these figures by giving an approximate B.C. (bolometric correction) 
for an 05 star. From Allen (1955).? BC = -42.5 + 10 log T + 29,000/T, and for 
T = 79,000 °K, BC = 6.9. Using the value, ^ = 4.62 for the sun and 

L^ 3°86xKr 'ergs/ sec. , we find 1.9xKr ergs/sec. This figure indicates 

that our calculated luminosity beyond the Lyman limit is of the correct order 

of magnitude. 



Within 11 parsecs, we have approximately 11 (05) stars, so that the total 
energy output absorbed by 2000 O of hydrogen is 1.9 x 10^® ergs/sec. Each atom 
of hydrogen absorbs 13.54 ev in order to ionize. In 2000 O, there are approx - 
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imately 2.4 x 10 hydrogen atoms* This means we requires 

2*4 x 10^ x 13 ° 54 ev x h6 x 10 ^ ergs/ev - 5-2 x 10^ ergs 

58 

to ionize all the material* At the energy output rate of 1*9 x lCr ergs/sec <, 
the ionization lifetime iss ^ qq 

b.z *10 e^c ,5 

|.q< IO^ g ero j s /sec. ^ 3 WO 7 

Such a short ionization time (this is., of course, a minimum limit because 
of neglect of radiation losses) would suggest that the interstellar material 
is almost immediately ionized as it appears and quickly -lost from the cluster. 




However, all of our estimates have been based on the use of Stromgren spheres 
for a group of stars with the interstellar material uniformly distributed. In 
reality, we would not expect this to be the case. Space densities derived for 
M3 from Sandage's counts for stars brighter than M.y= +3° 5 show that these 
brighter stars have a distribution proportional to R~^„ King (1962) has 
recently given an empirical law for projected, densities of the inner regions of 
globular clusters substantiating this type of distribution. Assuming that the 
initial space density of the material follows the same law as the brighter 
stars, we conclude that the initial distribution of matter would be highly con- 
centrated toward the center. If we also assume that the few stars responsible 
for most of the ionization follow the same R y law, then we can see that the 
central regions provide the likeliest environment for breeding HI clouds. The 

volume of ionization varies directly as the total luminosity of the ionizing 

p 

stars and inversely as N „ Since the total luminosity varies directly as the 
number of ionizing stars, we finally can write that the percentage of ionized 
volume varies as . In other words, the outer regions tend to ionize 

more easily. The foregoing arguments do not attempt to provide any mechanism 

to account for the actual formation of HI clouds, but merely offer a general 
analysis as to the most favorable regions for such formation. 
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An investigation by Spitzer (1942) shows that for a spherical self- 
gravitating system the interstellar material retained by a potential field would 
tend to concentrate toward the center because of the field and the effects 
of inelastic collisions. This effect would accentuate the initial concentra- 
tion. The fact that Roberts’s clouds are found some distance from the center 
seems to pose an anomaly. This point will be elaborated further in Chapter IV. 

It is evident from the foregoing discussion that one of the biggest 
detriments to the formation of HI regions is the presence of hot., blue 0 stars 
with their high proportion of radiation beyond the Lyman limit. This would 
suggest that a correlation exists between the presence and distribution of 
0 stars and the formation of HE regions. Any such correlation would not be 
straightforward. ,, however., because once dense HI clouds would form* then the 
subsequent appearance of these sub luminous 0 stars in the process of evolu- 
tion would not necessarily ionize dense regions. Nonetheless., it would seem 
logical that some correlation between the age of stage of evolution of a 
cluster and a favorable environment for the production of HE regions might 
exist. This point will be discussed further in Chapter V. 

We next look briefly at the kinetic equilibrium temperatures we might 

expect to find in the HI clouds of the cluster. A series of papers by Spitzer 

and Savedoff (1948., 1949^ 1950.9 1955) have concluded that the equilibrium 

temperature of the HI clouds In the galactic plane is nearly independent of 

6 

ins density in the range from about 1 to 10 H atoms cm . This conclusion 
assumes that the gas has sufficient time to reach equilibrium with its asso- 
ciated radiation .field. 

Before analyzing the HI regions in a globular cluster^ we first might 
discuss some of the more important energy gain and loss mechanisms that are 
involved. These were discussed by Eddington (1930) and later., in more detail 
by Spitzer and Savedoff. Throughout the discussion, we should remember that 
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we are dealing with a medium not in thermodynamic equilibrium, and in a dilute 
radiation field. However , the frequent elastic encounters between the particles 
do tend to establish an equipartition of the kinetic energy of translation, and 
the resulting approach to a Maxwellian velocity distribution can be defined, 
by an effective kinetic temperature. 

We can assume that there is virtually no radiation in HI regions beyond 

o 

the Lyman limit, 912 A. The primary heating mechanism is the superelastic 
collisions of electrons with those atoms having ionization potentials below 
13.54 volts and which are relatively abundant, i.e. C, Mg, S, Fe, Al, 

Heating by cosmic rays, at least in the solar neighborhood, is probably negli- 
gible for HI regions of 1 atom cm. or denser-. ’The energy loss mechanisms 
are (l) inelastic collisions of electrons with ions (particularly CII and Sill), 
(2) inelastic collisions between H atoms and grains (3) inelastic collisions 
between electrons and H2 molecules, and H atoms and H2 molecules. 

While Spitzer and Savedoff do not consider time “dependent sources, Kahn 

(1955/ estimates that collisions between separate clouds could be even , 

more effective than photoionization of the atoms as a. heating mechanism. By 

their analysis, Spitzer and Savedoff arrive at an equilibrium temperature 

o -3 

between 30 and 100 K for HI clouds of density 10 cm or more. At low 
densities between clouds where radiative cooling processes are ineffective, 
temperatures upward of 500 '"’k are postulated. Kahn shows that cloud temper - 
auures can be expected to drop very rapidly down to 500 °k between collisions, 
and more slowly below 500 °K. These analyses are Qbservationally substantiated 
by the 21-cm measurements which yield a harmonic mean temperature of about 
100 to 125 °K for the HI clouds of the galactic plane. 

In attempting to apply any such analysis to HI clouds in globular clusters, 
we are at once struck with the great uncertainties involved. Certainly, we can 
generally assume that the chemical composition of the medium in the clusters has 
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a lover proportion of the heavy elements than does gas in the galactic plane, 

Spectra of the globular clusters are interpreted to indicate metal to hydrogen 
ratios from l/lO to l/lOO that of the sum A lover percentage composition of 
the elements with low ionization potential would reduce the contribution of one 
heating mechanism for HI regions and perhaps lead to lower equilibrium tempera- 
tures. 

We can also expect the radiation density inside a globular cluster to differ 
from that of the solar neighborhood. We can estimate this value quite readily. 

We use the following notation in calculating the approximate radiation density 
inside Messier 3s 

u = radiation density ( ergs/eir^) (> 
r = distance from center of cluster ( cm) , 
c = velocity of light (3 x 10 10 cm/sec). 

T = average luminosity (ergs/sec) of those cluster stars brighter than +3° 5° 
L = total luminosity of cluster (ergs/sec) 

*/( r) = density of stars in cluster brighter than +3,5 ( stars/parsec^) 
at a distance r. 

From the Oort -van Herk data for model (d) of M3* we use the value 3,3 for L/qQ and 

190*000 solar mass for the total mass (0^) of the cluster. Using the value 

3-86x10^ ergs/sec* for the total luminosity of the sun and neglecting the probably 

small bolometric correction for the cluster* we find: 

L = (3-86 x 1 q 33 exgs/sec) ,( 190*000) , (3-3) = 2,43x10^ ergs/sec. 

For the purpose of determining the radiation density of the cluster* we consider 

that most of the light comes from those stars brighter than M^= +3,5* so that the 

density distribution of these stars can be equated to the distribution of light 

sources. For a single star: u = ^ at any point a distance r from the star, 

4jrr c 

Wow* to account for the contribution of radiation of all the other stars to 
this single point* we integrate over all the cluster* or in this case out to 
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100 pc,, the visible limit. 



So., at a single point in the cluster , we have 

lOO 



o) «= r: 



q-TTr^c 



fy-rrrW) dr - f 



loo 




For the total luminosity of the cluster,, we have the following expression; 

Urrr^V dr 



(?) ~ r? 



Dividing (l) by (2)., we get 



u 



J /.= 



j; 8 °vd. 



Wtt'c f o lo ° ir z V dr~ 



Our final expression for radiation density is then 

(1~ L ,€> ° Vd r 



(3) 



'47TC 



j^T^Vchr 

Figures (2) and (3) from Oort -van Herk depict the space densities (Vj of the 
group of stars brighter than M^= +3«5 versus cluster radius (r). These graphs 
are based on star counts for Messier 3 hy Sandage and von Zeipel. In order to 
find u,, we evaluate equation by carrying out the integrations. We split up the 

interval from 0 to 100 into 3 segments; (l) integrating numerically by Simpson's 

L 

Rule from 0 to 4 psc (2) represent analytically as a/ = — — 1 from 4 to 7 parsecs } 

, , , * 3 r 3 

and (3) as /y = — from 7 to 100 parsecs. 

Upon evaluating the two integrals , and using L = 2. 43x10^ ergs/sec,, we 
find that u = 2.9 x 10“^ ergs/em^. 

This value can be compared with the density of dilute black body radiation 

as calculated for the solar neighborhood; 5° 24 x 10”"*"^ ergs/cm^ (Dunham,, 1939) 

-13 3 

or 7 . 67 x 10“ ergs/cnr (Eddington,, 1930 ) } for which the black body temperature 
of a particle in space is 3«l8 °K. 
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4 -15 / 5 4 

Using the formula, u = akT , where a = 7» 57^10 ergs/ cm deg and k is the 

adsorption coefficient, we can solve for the effective temperature of the dust 
in the cloud if we can find an expression for k. 

We know from Roberts that the 260 O cloud absorbs approximately 5 magnitudes 

k , 29 

in the visible light ( X^xlO’^cm) . From Wien's displacement law, X " ' =1— • 

max l 



We can solve for 



X for the effective temperature, 
max 



We also know that k varies as l/x, so k = 2 . . . 



We therefore have k = " .,4,'- -..-ft 

.29 



r 

max 

25x10 5 

, giving us u = a T . Substituting 



our value of u = 2.9^10 ergs/ cm^ and solving for T, we find Tcs; 34 °K. 

The observational evidence supplied by Roberts (1959) by his 21 -cm measure- 
ments of M13 implies an equilibrium temperature of 100 °K or below. Hence, we 
are led to conclude that the temperatures of dense HI clouds in globular clusters 
are at least as cool as HI clouds in the galactic plane. 
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CHAPTER IV 



PASSAGE OF GLOBULAR CLUSTERS THROUGH THE GALACTIC DISK 

We now turn to an investigation of the interaction Between gas clouds and 
glohules in the cluster with the medium lying in the galactic plane, Previous 
investigators have simply stated that any gas present in the cluster will he 
ejected hy each passage through the disk. 

The galactic disk consists of a thin layer of gas and dust distributed in 
a non-uniform fashion. There are individual clouds and spiral arms separated 
hy regions of low density. Densities for each feature ares typical individual 
clouds - 10 atoms cm ^ spiral arms - 1 atom cm * space between clouds - 0.1 to 
0.01 atoms cm , As previously argued* we assume most of the hydrogen to he 
neutral for purposes of the collision process. The ionized gas in the disk is 
concentrated around the individual 0 and B stars. 

Since we are interested In the processes that tend to act on the neutral 
hydrogen clouds and glohules In the clusters* we can see that the collision to 
he studied is the HI -HI impact. 

As long as the mean free path for particle collision is small compared to 
the dimensions of the cloud or globule* we can treat the gross dynamics of the 
encounter hy evaluating the momentum transfer to determine if the clouds and 
glohules of the cluster are ejected from inertial force considerations. For 

o 

collisions of neutral H atoms* we can use the geometrical radius (~10“°cm) to 

if) P 

compute the collision cross-section, giving CF= 3x10“ cm (Parker, 1958 )- For 

r 3 

the globules, N — Krcm % for the clouds, N ‘^lO'cm . Using the relation, 
rL n 

M mean free path)= l/N<T, we find \ glob= 10' 8 pc and \ cloud— 10"^pcj with both 
values considerably smaller than the dimensions of the object involved,. 

Kahn (1955) has considered the statistics of size and distribution of clouds 
in the galactic plane by using the observational evidence that: (l) a line of 

sight of 1 kiloparsec intersects an average of 10 clouds, and (2) clouds occupy 



- 21 - 



about 5$ of the volume of the galactic disk. He concludes that the average 

cloud is 3*75 parsecs radius, with a density of 20 H atoms cm and that the 

J 4 -3 

number density of clouds in space is 2.45x10 pe „ This result, of course * is 
characteristic of the solar neighborhood. We readily see that the mean free 
path for collision is about 100 parsecs. 

From 21 cm. observations, Rougoor and Oort (i960) have described the general 
galactic disk between 5 an( l 10 ^pc from the galactic center as being 220 pc thick 
with a mean density of 1 H atom cm The disk inside 3 kpc seems to be much 
thinner, about 120 pc, and the mean density is given as 0.4 H atoms cm“^ 0 They 
comment that the region between 3 kpo and the central disk at 600 pc appears to 
be virtually empty. 

From our discussion in Chapters II and III, we note the following informa- 
tion concerning the clouds and glob -dess (l.) In general, we would expect to 
find HI regions toward the center of the cluster in an isolated cluster. 

(2) From the incomplete information now available, clouds seem to appear at the 
edge of the clusters (minimum of 8 to 10 pc) and globules toward the center of 
a cluster (in the one example, M4). 

Let us assume that a globular cluster contains a 260 O cloud of Roberts ? s 
description at the center of the cluster. Let us then analyze the possible 
effects of a collision with the galactic disk. In this analysis we will use 
two different treatments for the two different types of impact involved. If we 
assume that the cluster traverses a mean distance through the disk of 300 parsecs, 
then at 300 km/sec (where 1 km/sec C^l parsec/l0^ years) we can conclude the 
time of transit is about 10^ years. However, a cluster cloud would transit one 
of Kahn’s typical, dense clouds in the galactic plane (diameter^ 7° 5 P c ) in 
about 10 r years. If we compare these two characteristic times with the approxi- 
mate orbital period of an object in the potential field of the cluster, 10 years, 
we see the need for different treatments. 
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Momentum Transfer 



For the case of impact with a discrete cloud in the disk, we consider an 
impulse to he acting on the cluster cloud- We will also assume that the colli- 
sion is completely inelastic, the colliding portion of the disk cloud remaining 
with the cluster cloud. 

For the case of impact with the intercloud medium (0.1 to 0.01 H atoms cir ‘" > ) 
of the disk, we compute the momentum transfer per unit time of the colliding 
medium, which is the force acting on the cloud or globule of the cluster. 

We will see that one effect of the collision with h disk cloud can he to 
knack a cluster cloud out of the cluster center and into a rectilinear orbit. 

On the other hand, the force exerted by the medium on a cluster cloud effectively 
shifts the whole orbit in the direction away from the direction of motion. 

The simplest manner in which to visualize this problem is to analyze it 
with respect to the potential well of the cluster and to consider the motion 
as a one dimensional oscillator. We use the following notation; 

V = velocity of globular cluster relative to the Galaxy (300 km/sec). 

R = distance of cloud or globule from center of cluster. 

r = radius of cloud or globule of cluster. 

L = diameter of disk cloud ( 7*5 P c ) 

n = number density of H atoms in galactic disk (cm'”^). 

-24 

m = mass of hydrogen atom (1.67x10 gm). 
ri 

(V) = mass of cloud in cluster (260 O). 

<$(R) = potential as a function of R in the cluster (km/sec) . 

= potential gradient. 

N = number density of H atoms in globular cluster cloud. 

F = Momentum per unit time or force applied to cloud in cluster resulting 

from transit of intercloud medium of disk. 

E = Total energy per unit mass of cloud in globular cluster. (km/sec )^. (with 
respect to cluster frame of reference). 
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Figure (l) depicts our detailed example of the globular cluster with a 260 O 
cloud at the center of the cluster on its first impact with the disk. 

First, we derive a general formula for computing the force per unit mass 
resulting from momentum transfer per unit time from the intercloud medium to 
the cluster cloud on transit of the disk. We have that: F = nr .V' . n.m_,( dynes) 

3 ** 

09 = . N. m H ( grams ) 

(1) F/fl7= (cm/sec 2 ) 

For our- particular case, where r = 1 pc and N = 2.7»10^cm 1 V = 300 km/sec and 
n = 0.1 cm~f equation (l) gives us: F/fl? = 8. 2x10*"^ cm/sec 2 or 2.5^~fc^ _ P- t 
We can also give for our example, the expression for the force as: 

(2) F( dynes) = 4.25xl0 23 . n( cm" 3 ) 

or F = 4.25xl0 2 ^ dynes (for n = 0. l). 



We have the following relationships between total energy per unit mass (E), work 

F.(E -R ) 

done on the cluster cloud per unit mass ■ , — potential energy per unit 

(u 

i 

mass ($) and kinetic energy per unit mass (5^), whereof - velocity of cloud in 



its orbit with respect to the globular clusters 

*p ^ 

(3) E = E_ ) where E is the total energy in an orbit at time t^ 

d 1 Vi ^ 1 1 

and E^ is the total energy at some time t^ later. As shown in figures (l) and (2)^ 
R values in the direction away from the direction of impact are positive.. We 
also have the following identity relating our energies s 

_ i 2 

E - (where HT is the velocity of the cloud in its orbit relative to 

the frame of reference of the cluster) 

V 2 

We can see that the value of f/m «= —4 determines to what extent an orbit is 
displaced to the rear because of the steady force exerted by the intercloud 
medium during about 10° years . In figure (2) we have depicted the effect of 

~z o i 

force per unit mass exerted by both a 0.1 ciri medium (f/ 0Q = 2.5(km/sec) pc” ) 
and a 0.01 cm ^ medium (f/|Y| = 0.25)* We depict the extreme case for both these 



- 24 - 



media in order for a cluster cloud (cf Roberts's description) to be retained 
in the cluster for an indefinite period of time. For the 0.1 cm medium., the 
limiting orbit is from -4 pc to +14 pc; for n - 0.01 cm“i it is from -20 pc 
to +50 pc. 

We should note , however, that for oar particular case where the time of 

6 7 

transit of the disk is 10 years and the orbital period in the cluster is ~10 

years, that the chance of a cloud escaping due to the reaction of the medium alone 

is extremely remote,* Once the cluster emerges from its transit of the disk<> 

then the f / Yj slope is practically zero and^ as long as the cloud has a total 

o 

energy less than 138 (km/sec) , it will not escape. 

However,, the reaction caused by collisions with discrete clouds in the 
galactic disk of the order of n = 10 cm ^ is extremely significant for our 
discussion* To calculate the energy change per unit mass (AE)^ we use the 
principle of conservation of momentum. For the general case we have for the 
mass of material intercepted by the cluster clouds 

Lo jt. n.m^ 

The mass (OT of the cluster cloud iss 



V is the relative velocity of impact (300 km/sec). 

AHTps the velocity change imparted to the cluster cloud relative to the cluster 
frame of reference. 

We then have; PQ;.V = (09 +P9,-).A'*T 
This gives the expression ATT c 



(4) A0T= 



V 

l"+ 4rW 

3Ln 



Now we derive a general expression for the energy transfer per unit mass per 
collision. We have; 
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i 2 

E_^ = $+ 2^1 (where f \T_^ is the orbital velocity of the cloud at time t^) 

E = $+y( d/7 +A (\f) r ~ (E^ = energy after impulse) 

This gives us the general expression for energy transfers 

2 

(5) AE = E -E ± = nr AV+ L^L 

Substituting the appropriate values into equation (k) for the 260 O cluster 
cloud colliding with the n = 10 cm"^ disk cloud, we find A die: 6 km/sec. From 
equation (5) then, we have AE = 6 (dT. +3 ) . 

Wow, with the use of figure (l), let us follow through in some detail the 
reactions of a 260 O cloud, initially at rest in the center of a globular cluster 

as it passes through a typical area of the galactic disk. We adopt 0.1 H 
-3 

atoms cm ■ as the concentration for the intercloud medium. We assume that the 
average cluster passes through .300 pc of the disk. The mean free path for 

collision with one of the disk clouds (L = 7° 5 n = 10 cm~^) is 100 pc. 

~3 

Upon gradual entry of the cluster into the 0.1 cm v intercloud medium, we 
expect the cloud to be shifted about 0.2 pc to the rear of the cluster center, 
or to R = +0.2 pc. The cluster cloud is likely to be near this point on the 
diagram upon its first collision with a disk cloud. Using our previously calcu- 
lated value of AU^ 6 km/sec, we find from equation ( 5 ) , where'll^ 0, AE = l£ 

We indicate this change in energy on figure (l) by a dotted vertical line to 
point (R^, E ) or (+0.2, 19)- At this point we assume that about one third of 
the transit period ( ^xlO^years ) has elapsed. Through this point we draw a line 
with the proper slope to represent F/fQC£2,5 (~ c ) 2 pc“ 1 , as previously calculated 
for the force per unit mass exerted by the 0.1 cm“^ medium on the 260 O cloud. 
This line represents the locus of (R,E) points for the rectilinear orbit which 
the cluster cloud now follows. 
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The increment of energy., AE = 18 (-8—) , received by the cluster in the 

S 6 C 

first collision has given the cloud an impetus toward positive R values * We 

now estimate the distance traversed along this orbit in a time interval* 

5 

3x10^ years* until the next collision* We integrate numerically along the 
path from R = 0« 2 pc tc an R value that corresponds to a time interval of 
3 x 1 Ck years. We use distance steps of 0.2 pc; 

cffl 

\ e x ?0TeT 

We then find that the cloud will he at (+1.6. 22) at the time of second collision. 

As depicted in figures (l) and (2), we follow the cloud through two more 

s 

collisions. After the second collision and a time interval of 3x10 years, 

'2 

we find that a AE ,= 58 (ton/ sec) for the third collision raises the cloud to 
a very high energy total (e ss 125 (ton/see)^ which means the cloud probably 
escapes. 

We note from figure (l) that if the cloud does not undergo the third 
collision, it will find itself in an orbit, with the distance limited to about 
7 parsecs. From figure (2), we see that a cloud after the third collision is 
supposedly left in an orbit extending out to about 60 to 70 parsecs. However, 
it is easy to see that, even for a very small frictional force which could 

-2 -4 “3 

reasonably be expected from the halo medium (10~ to 10”' cm ' ), the cloud 

would be lost from such a high energy orbit. 

While we have worked out only this one example in any detail, it is easy 

to judge the possible reactions in other situations. 

For instance, a subsequent pass by our cluster in the example above^ which 

8 

has had two collisions on the previous pass through the disk 10 years before ^ 
might find its cloud in an orbit at an energy level of about E = JO* The 
initial velocity* for this case might vary anywhere from about = +12 km/ sec 
through 0 and back to -12. Minus values for velocity indicate a direction for 
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..ue cloud in the same direction as the motion of the cluster orbit. Using 
formulae (4) and ( 5 ) , we find that AE can have the following values for each 
value of V 1 for +12 km/sec, AE = +90? 0 km/see, Ae = +l8j -12 km/sec, Ae = -54, 
Hence, we see that in these 3 cases," the cloud escapes, has its energy raised 
about 18 (km/sec) 9 and has its energy lowered 5^ (km/sec) respectively • 
Statistically* we would expect the cloud to be near an extreme point of its 
orbit and thus find that the change in energy for an impulse is small 0 

Two other factors are worth brief mentions (l) we,>could expect some 
elliptic orbits to ensue^for the reason that a cluster on a subsequent pass 
could collide with the disk at an angle to the direction of motion of its rectilinear 
orbit j and (2) we might expect some damping of the orbit during the period the 
cluster is away from the disk* because of accumulation of additional medium 
from evolution of stars and the effects of inelastic collisions 0 

To estimate the chance of a particular cluster cloud getting through the 
disk without escaping from the cluster* we use the results from the previous 

discussion and apply the Poisson frequency function: 

k -a 

f(k) « 

where a is the expectation and f(k) is the probability that x = k. For our 
problem, we take a equal to 3 (we can expect about 3 collisions per pass). 

We then evaluate the individual probabilities for k = 0, 1, 2 and 3 collisions 
per transit. We sum up these individual probabilities, using f(k) for k = 3> 
and find that the probability of a cluster cloud surviving one disk passage is 
about kc Of course this figure is for the average assumed conditions of the 
disk and would vary for a different structure in any region. 
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We should note that the foregoing development for momentum transfer 

applies equally well for the Idlis cloud. That this is true is readily 

understood from equations (l) and (4) which show that the forces depend 

on the product rW which are approximately equal for the two cases. 

It is clear that any original gases in the globular cluster are 

10 

swept out at some time during passages through the disk during 10 years. 

We previously mentioned one other possible source of material for the cluster., 
that of sweeping up the ga.se s from the halo. Spitzer gives a value of 
5x10 “ h auoms cm“ for the halo. Shklovsky and Pikelner (1958) give 

-7 

10 cm . If, for example, a Roberts cloud swept a path during its orbit 
of approximately 30 kpc length, and retained all swept matter. The total 
mass accumulated would, be only about 20 O, even by using N = 10 We 
can conclude that a "snow plow" accumulation is a relatively small source 
of material. 



- 31 - 



iergy Transfer 



In addition to the effects of momentum transfer^ ¥e must also consider the 
heating effects of a collision between clouds at these high energies* From our 
discussion in Chapter III of heating and cooling processes in cosmic clouds., we 
learned that if an HI region becomes ionized^ it will have an equilibrium temper 
ature of about 10^000 °k and the root mean square velocity of the protons will 
be about 1 6 km/sec* The Oort-van Herk dynamical model for M3 gives a velocity 
of escape of 1 6*6 km/sec at the center ranging down to 11*4 km/sec at 8 pc and 
so on^ as can be surmised from figure 2. We conclude that when hydrogen is 
raised to a kinetic temperature of TO .>000 it will quickly expand and escape 
from the cluster* 

Bates and Griffing (1953) have calculated the ionization cross section for 

the reaction H(lS)+H(lS) = H(lS)+iI + -he, Extrapolating their data^ we estimate 

20 2 

this cross section as 10 cm for a velocity of about 300 km/sec* Comparing 

oo o 

the ionization cross section (~10“ w cm ) with the collision cross section 
/ 16 2\ 

V~lu cm ') , we see that an individual hydrogen atom has a very low probability 

\ 10 ) of ionization per collision. The incoming energetic hydrogen atoms, 

relative to the frame of reference of the dense (~10^ to lO'f rnU ) cloud or 
globule of the cluster, will be quickly damped to energies below 13.6 ev as 
a result of both momentum transfer and radiation losses by excited atoms. As 
we will see later, the ionization will he mostly confined to a narrow interface 
between two colliding clouds. 

We should keep in mind that collisional ionization is not a very efficient 
process. Alfven (i960) has proposed that collisional ionization between a 
plasma and a neutral gas is caused by electrons being accelerated to sufficient 
energies by irregularities in the charge distribution created in the plasma. 

In the discussion following this paper, Spitzer remarks that when ionization is 
produced by energetic electrons, the cross section for excitation is much greater 
than that for ionization and h to 5 times as much energy goes into light than 
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into ionization. Naturally., for any ionization mechanism* a certain proportion 
of the input energy will be dissipated as light and infra-red radiation from 
excited atoms and from grains. 

We will analyze in some detail the case of collision between a Roberts 
cloud (N~10^em”^) and a typical cloud in the disk (n~10 cm”t diameter [l]~7„5P )• 

First* we should note that at relative velocities of about 300 km/sec, the 
interacting particles exceed the velocity of sound and we can expect a compressed, 
shock front between two colliding clouds. The velocity of sound is given by 
a =\[^ RT* where 5/3* R = 8. 3x10 t For T ^ 100 ’“K* a = 1.2 km/ssej for 
T = 10*000 °K* a = 17 km/ sec. 



We now calculate the total energy per second transferred to the more dense 
Roberts cloud as it passes through the larger ‘loud of the disk. We have for 
energy/unit time; 



E/unit time = Jtr^xV )x( rn . nrg).xV^ ~ 6. 4x1 0* trgs/sec 

(where = 10 cm“^ for disk cloud* r = 1 pe for Roberts cloud* V = 300 km/sec* 
'W 1.67x10“-" gm). This means that ; or 0.23 ergs cm~ sec is the 



jtr 

amount of energy input per unit time per unit area which must be dissipated at 



the interface of the Roberts cloud. We note that this value is comparable to 

the average stellar radiation per unit area in a cluster \ 

uc 2.9xl0~ 10 x3xl0 10 . -2 -1 

= — - — ^ ■ — - — = 0.22 ergs cm sec 

We should investigate the possible cooling mechanisms which are able to 
dissipate this energy and maintain the layer at a reasonable equilibrium temper- 
ature , 



Spitzer ( 1949 ) has given the following formula for computing the rate of 
energy loss caused by the inelastic collisions between hydrogen atoms and dust 

grains; 

L ( ergs/cmtee) = 9- 4x10 rii g)N(H)T :5 (T-Tg) 
hg 
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We take tne following values tors 



O -z 

n(g) - concentration of dust particles = 10 “cm~^ 

C Z 

(where radius of particles Qfl.O“ v cm and density csl gm/cnr ) 

W(H) - value for Roberts cloud^lCr cm~^ 

(the above two values give a mass ratio of dust to gas of approximately 1 fc' 
T - gas temperature assumed to be 10., 000 '“K. 
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T - temperature of the grains is assumed very small (10-50 K ) . 

§ 

01 

Inserting these values in Spitzer's formula,, we find L = 9» ~xl3"' ergs cm ^seo’ ", 

Hg 

We have an energy input by collision of 0.2.5 ergs cm~2 seo ,“'^ which even if dis- 
tributed throughout a layer of several mean free paths fox collision 
(1 = l/w<T = 1/lvl ,5x10 - 3xl0 12 cm) would be many orders of magnitude greater 



than the cooling rate. 

We assume that- the impact- layer between the two clouds is completely ionized 
and at a temperature of 10,000 Accordingly, we will have a recombination 

of electrons and protons with the accompanying radiation in the visible and. 
infra-red as the electrons cascade down to the ground state. We use the fellow- 



ing formula to obtain an estimate of the cooling rate provided by this mechanisms 

Qo) L (ergs cm ~see~"*") = N N /IT 'y-tT ( h •s/ ) .D 
rec® p e n n 

z 

where N = proton concentration in the impact layer (cm ) . 

N = electron concentration (cm )„ 
e 

Of" = velocity of electron at 10,000 °K® 

(f m = capture cross section to n^ 1 state for hydrogen® 

= average energy radiated per capture to state ri, and in cascade 
to ground state® 

D - thickness of ionized impact layer (cm)® 

From page 9 1 Allen (1955)* we take the values given for the recombina- 

th 

tion cross section to the n ' state for hydrogen and multiply each value by 

(hV) r and then sum all products to give; ^(T^hV) = 5. 7x10 "^ergs cm”^ 

(where each 0*^ is of the order of 10^~cmt . 

= 3k - 



For the velocity of the electron at 10 ,,000 °K, we 

consider that N — u and set the equation for L 

p e rec 

We then solve for N 2 ^ 



use (IT= 6. 4x10^ cm/sec. We 
equal to the energy input,, 



N?D(V)(^fT (h-v/) n = energy input from collision 

ift) (6.4x10' )(5. 7X10 ~ 32 ) = 0.23 ergs em^sec' 1 
( 7 ) or ITD = 6.3xl0 22 cm" 5 

O 

Having obtained an expression for N we now find N from considerations 
of pressure equilibrium in the collision process * We note that the collision 



picture can be described in the frame of reference in which the ionized shock 

layer between the two clouds is stationary so that a disk cloud (n = 10 cm~^ ) 

moves with a velocity of about v = 300 km/seco We can readily compute the rate 

of transfer of momentum per unit area at the front by the expressions 

n o m^ 0 v^- 10 ( cm~^) • !• 67. 10^"*( gm)o (3xl0^)^cm“see"^ 

i r vr 8 2 ~ 2 

- lopolO gm cm sec 

Since this value is the force per unit area or pressure exerted on the 

layer from the direction opposed to the direction of motion., we can equate it 

to our expression for the pressure inside the layer^ 

P = K kT (where N = number of protons per cirT% k = Boltzmann’s constant 
P p 

(1,38x10“ ergs/degree) ^ and T — lCpOOO °K). We use N as the concentration in 

P 

this expression since the electrons contribute only about l/43 of the effective 

II 

pressure. We find N ctlO cm C 

P 

Using our result from equation ( 7 ) that N^D = 6.3xl0 22 cm*h we then solve 
for the thickness of the shock layers finding D = 6. 3xl0"*"^cm. We note that D 

"I Q 

is larger than the mean free path for collision of H atoms ( X = 3x10 cm) in 
the dense cloud and comparable to the mean free path (X =3xl0^ ,+ cm) in the disk 
cloud. We also note that it is extremely thin compared to the diameter of the 
cluster cloud (2 pe). 
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It follows that we have a situation in which the two clouds and their 



interface are in pressure equilibrium. Using the expression for momentum 



2 P 

transfer per unit area., n.ir ,v , we see than s Or . 

h ’ u cluster 



V 2 n (300) 2 .10 
W = 10 3 



or ^cluster 



30 km/sec (relative to frame of reference where shock layer 
is stationary). 

We then see that the distance of travel by the shock front relative to 
the cluster cloud during the time of collision is equal to -^2, x L (where 
L = 7-5 pc). This distance of travel by the shock front ("-0.75 pc) is less 
than the diameter of the cluster cloud (2.0 pc), as we would expect. 

The effect of this impact- layer is to convert high energy particles int-c 



radiation in the visible and infra-red region at the spectrum.. 'This., in turn., 
serves to increase the radiation density at these wavelengths inside the cloud 
of the cluster. In order to determine if the cloud can 'handle this increase in 
radiation density at these wavelengths, we analyze the cooling capacity of just 
one mechanism - that of the dust grains. If we should, find that an increase in 
radiation density would require such a high temperature for the dust particles 
that they could not remain in the solid state, then our cloud could possibly 



be raised to a high equilibrium temperature. We pursue this analysis in a way 

which seeks merely to set a limit as to the amount of cooling the dust can handle 

with respect to the increased energy input. 

From Roberts 1 s data we have that the cloud absorbs about 5 magnitudes in 

the visible region of the spectrum (X~3xl0" 3 cm). If we assume that the dust 

has a temperature of about 50 °K (at 50 °K, the cloud will be self -gravitating), 

then from Wien's law: X. (cm) = — 

max. T(°K) 



X (at 50 °K ) ~ 5 . 8x10 ” 3 cm 
max x ' 

Since the absorption law for small (~10“5 em ) dust particles shows that absorp- 
tion varies inversely as wavelength, we conclude that- effective absorption 
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0.05 magnitude. This means that 




at 50 K is approximately 5 mag| 
the cloud has an absorption coefficient, (k)s£5$, at the relevant wave- 



length , X 



max. 



Therefore , we apply the Stefan -Boltzmann radiation law (with the additional 
factor k) to our cloud, and thereby determine hew much energy per second the 



dust in the cloud can radiate at 50 °K; 

L(ergs/sec) = 4xtr c '<r T^.k 
where r = 1 pc, <T = 5.'7xlO"^erg cm“ 2 deg “^see" 1 , 

T -.= 50 C K and k = C.05. 

Substituting these 'values, we find that L~ 2xlC’ ergs/sec ( coding rate that 
can be handled by a Roberts cloud, with dust particles at- a temperature of 50 °k). 
Now if we compare this with the maximum possible energy input from the collision 
process, 0.4x10' " ergs/ see, we conclude that, the cloud can absorb this addi- 
tional energy input without vaporizing the dust:. 

By comparing with the aboveythe energy input expected in the collision of 
a luster cloud with the intercloud, medium of 0.1 cm v concentration, we con- 
clude that the Roberts cloud will not- be heated, appreciably except for a thin 
layer at the collision juncture. 



We also find that the Idlis cloud can withstand the thermal effects of 
cloud impact in the galactic disk. 

Ey a similar analysis as used for the 260 O cloud, we find that the 

k ^ lb 

proton concentration is also 10 cm in the Idlis cloud, making D nr 6x10 cm. 

Using the expression for transfer of momentum per unit area, we find that the 

ionized layer travels at 0.3 km/sec through the globule (of — - i ^ 22 ' ° 12 

globule q 5 

Since this is below the speed of sound (a = 1.2 km/ sec for the material} , 
we do not have a shock phenomenon in this case. 
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CHAPTER V 



DISCUSSION OF RESULTS AND CONCLUSIONS 

From the calculations in Chapter IV, we can conclude that a globular cluster 
will not necessarily lose its interstellar material in passing through the galactic 
disk. 

We have seen that, for a cloud of the proportions described by Roberts' and 
also for the smaller, denser globule of Idlis, momentum transfer will be sufficient 
in about half the cases to eject it from the cluster, We also learned that a 
shock front, forms doming collision with clouds of the galactic disk and. that the 
recombination of hydrogen^ and radiation in the infra-red by dust grains is a suf- 
ficient cooling process to dissipate the energy of collision which, therefore, 
will not cause heating and expansion of the cloud. 

From considerations discussed in Chapter III, we concluded that the center 
of the cluster was the most satisfactory environment for HI regions. 

Accordingly, it would seem plausible that diffuse material accumulates in 
f j center of the cluster and is excited into orbits at various distances on 
passages through the galactic disk. We would not expect to see those clouds in 
the cluster at great distances from the center, even if they occurred there, nor 
those near the center on photographs with exposures typical of those used by 
Roberts, 

We have examined at some length possible correlations for those 12 clusters 
( which show obscuring regions) out of the 32 examined by Roberts, Correlations 
could be expected for at least two reasons; (l) differences in interstellar 
avironments and evolutionary activity among clusters, and (2) significant, 
ystematic differences in the structure of the galactic disk. We must use caution 
u attempting a correlation since it is not at all certain that the 20 clusters 
without reported obscuration contain less diffuse material, Roberts has stated 
(private communication, 19 6 . 3 ) that he was conservative in his estimates of 
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clouds and that his list of 12 could possibly be increased by 5 0$* 

We have noted that only about half of the globular clusters have been re- 
ported as having RR Lyrae variable stars (Hogg* 1959 )» Adding co Cent aura 
(Fitzgerald 1955) and M4 (idlis) to Roberts* s list of 12 clusters* we note 
that RR Lyrae stars have been sighted in all but two* Of the 20 clusters 
examined by Roberts with negative results for obscuring clouds* 13 are listed 
u.s having RR Lyrae stars* Generally speaking* the lack of RR Lyrae stars in a 
cluster is attributed to the absence of stars In the app.rc-pri.ate evolutionary 
stage* While the above figures don’t indicate a, significant difference* in 
this ca.se* some similar criterion having evolutionary significance might possibly 
u uncovered by a closer search and classification of variable stars and extremely 
blue stars in globular clusters* 

We learned in Chapter IV that collisions with the dense clouds of the disk 
is the major factor- for lose of clouds fi xe the globular clusters* The force 
xerted by the intercloud medium is net significant* Accordingly* we might 
expect that if there exists a large* contiguous area of the galactic plane with 
a significantly smaller proportion ol dense (10 cm ') structure than in the 
solar neighborhood* that those clusters passing through this region in their 
galactic orbits would retain more interstellar material and be more likely to 
show obscurations* 

Bougoor and Oort (i960) have remarked that there appears to be a relatively 
small amount of dense HI structure in the region 600 pc to 3 kpc distance from 
the center in the galactic plane* Optical observations of the Sb galaxy 
:IGC 2841 (Sandage* 1961) show an amorphous region devoid of dust or spiral 
1 mature in the same region (as scaled to our Galaxy)* 

We have determined* by the use of a formula from Eggen* Lynden-Eell* and 
Sandage (1962)* that the average orbit of a globular cluster (e = 0*8) with 
apogalacticon, 10 kpe 5 sweeps out about 205° between successive apogalacticon. 
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"We also have determined that our average 0,8 e orbit will pass within 3 parsecs 
of the center of the Galaxy during disk transit if the semi -major axis lies 
within a cone with vertex angle of 50° ^ centered on the galactic axis. We 
have attempted a correlation, without success, between the present location 
of clusters within the cone and the presence of obscuring material. The fact 
that heliocentric distances of‘ globular clusters have probable errors up to 
.cakes any such correlation doubtful. 

Perhaps the two most significant recommendations for the future study of 
globular clusters are? (l) An observational program that will determine better 
the number of hot, blue stars in globular clusters. For example, the incam - 
■ lateness factor of about 15 for >5 stars in M3 could be improved by more care- 
ful observations. These stars are significant in their effect on the ionization 
of HI c-louds, and the radiation field of these stars might be the most important 
f actor influencing the life mime of HI clouds, since they have a reasonable chance 
•f surviving passage through the galactic disk. (2) It is desirable to have 
a more detailed program for observing its cured regions in globular clusters. 

Of primary interest would be the photographing of clusters at many different 
exposures to determine better the cloud locations and their absorption. 
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